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ABSTRACT 

We present the results of colour-selection of candidate high redshift (2.6 < z < 3.9) 
galaxies within the Hubble Deep Field based on the Ly-break at 912A. We find 8 such 
objects in the region, giving a comoving number density comparable to that of nearby 
bright galaxies (for a flat qo=0.5, Ho=100 kms^^Mpc^^ universe). We provide basic 
data on the properties of these objects, and show that despite their absolute magnitude 
being significantly brighter than L* (typically Mb = -22), they are generally smaller 
than nearby galaxies. Furthermore, visual inspection of their images shows that they are 
all highly disturbed systems, with multiple nuclei, tails and plumes, suggesting that they 
are undergoing merging processes similar to most nearby starburst galaxies. Theoretical 
models suggest that galaxies form by accumulation of numerous subcomponents, and 
we suggest that we are seeing this process underway in these objects. It is thus possible 
that the epoch of galaxy formation might have been discovered. 

Key words: galaxies: formation - galaxies: primeval - galaxies:interactions - galax- 
ies:starburst 



1 INTRODUCTION 



The discovery of very high redshift and primeval galax- 



tion, faint imaging surveys, for both number counts (Tyson 



1986) and to find individual interesting objects selected by 



various colour criteria (e.g. multicolour searches for high red- 



ies has been a 'holy-grail' for cosmologists for many years shift quasars) are well established techniques ( pmith et al 



( [Partridge 19741). Amongst numerous approaches to finding 1994! ). Among the more interesting results to emerge from 
forming galaxies and understanding their subsequent evolu- these studies has been the discovery, at blue wavelengths. 
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of many more galaxies than standard 'no-evolution' models 



predict - the so-called 'faint blue excess' (Tyson 1988). How- 



ever, various spectroscopic (Glazebrook et al. (1995a), Lilly 
et al. ( 1995| )) and HST-based morphological studies (Driver 



et al. (1995), Glazebrook et al. (1995b)) have shown this 



excess population to reside at only intermediate redshifts 
{z < 1), and to consist of dwarf late-type spiral/irregular 
galaxies which undergo substantial luminosity and colour 
evolution between z ~ 1 and the present time. On the other 
hand, these same studies have shown the redder early-type 
(E/SOabc) population to have undergone little if any evolu- 
tion in this period. The formation epoch of these galaxies, 
therefore, must lie at considerably greater redshift, even if 
the nature of the faint blue excess has now been determined. 

In this context, the Director of the Hubble Space Tele- 
scope (HST) assigned a substantial allocation of discre- 
tionary time (~ 150 orbits) to a project aimed at obtaining 
very deep Wide-Field & Planetary Camera 2 (WFPC-2) im- 
ages of a single 'random' field (position 12h 36m 49.4s -|-62d 



12' 58", J2000 (Ferguson et al. 1995); hereafter referred to 



as the 'Hubble Deep Field' (HDF)). Four filters were chosen 
for these observations, F300W, F450W, F606W and F814W. 
The data, which were taken over the 13-day period 1995 De- 



cember 18-30, were made public immediately ( Williams et 
al. 199£|) and we present here the results of a programme to 
search for high redshift counterparts of normal galaxies us- 
ing the presence of the Lyman break at 912A in the F300W 
filter as an indicator of redshift. 



discuss the properties of the objects we find. We draw our 
conclusions in section 5. 



2 SELECTION OF HIGH REDSHIFT OBJECTS 

The HDF extends to magnitude limits (on the STMAG sys- 
temQ) of 26.7, 28.8, 30.1, and 30.3 in the F300W, F450W, 
F606W, and F814W filters, respectively, for a 5 a detec- 



tion in 16 connected 0.04" pixels (Giavalisco et al. 1995). 
For the present study we wish to detect and then examine 
the morphological properties of high redshift galaxies. We 
identify high redshift candidates by looking for a substan- 
tial magnitude difi'erence between the F300W and F450W 
filters. For galaxies with 2.6 <z<3.9, the F300W passband 
will correspond to rest-frame wavelengths shorter than the 
Lyman cutoflt at 912A and the F450W filter to Ion ger wave- 
lengths. These high redshift galaxies will thus be unde- 
tectable at F300W or will have anomalously red F300W - 
F450W colours. Objects at still higher redshift will addition- 
ally have F450W suppressed. 

A similar method has been used successfully by Steidel 



et al (1992) to look for companions to high redshift quasars 



and damped absorption-line systems, and by Guhathakurta 



et al (1990) to set a limit on the redshift of the faint-blue- 
galaxy population. Additionally, the spectral energy distri- 
bution of objects at high redshift is likely to be fairly flat 



redward of the Lyman cutoff ( Djorgovski 1992 ), an assump- 
tion that is borne out by recent Keck spectroscopy of z~ 3 



galaxies (Bteidel et al. 1996). Our selection technique is thus 



In the next section we describe our selection technique. 
In section 3 we present our results, while in section 4 we 



* STMAG = -2.51ogio/A - 21.10 where fx is the flux in 
ergscm"^ received in the bandpass 



to look for all those objects with F300W - F450W > 2, and 
remaining colours (F450W - F606W and F606W - F814W) 
between -0.6 and 0.6. Consideration of the colour-colour di- 



agram in (Steidel et al. 1995) and the detailed photometric 



modeling of Fukugita et al ( 1995 ) demonstrate that it is very 



unlikely for lower redshift galaxies of any type to have such 
colours. Objects not detected in the F300W filter are treated 
as if they have a magnitude of 27.7, corresponding to the 2a 
upper limit, for the purposes of this selection. We are thus 
constrained to objects with an F450W magnitude of 25.7 or 
brighter for our candidates, which usefully guarantees that 
they will be bright enough for us to study their morphology. 
Examining the flux zero points for these filters shows that 
for an F300W - F450W colour > 2, the fiux ratio between 
the F450W and F300W bands is > 6. This is higher than 
could be caused by the 4000 Abreak or by absorption in the 
Lya forest where the continuum suppression factor at these 
redshifts is only about 50% ( ^torrie-Lombardi 1994 ). 



This selection was based on object catalogs derived from 



the combined 'drizzled' HDF images ( Williams et al. 199(: ) 
using the automated photometry program SeXtractor writ- 
ten by E. Berlin. A detailed description of the techniques 
this program employs in the detection and photometry of 



objects is given by (Berlin 1996); its specific application to 



HST WFPC-2 images is discussed by (Small et al. 1996) 



In using it here, SeXtractor was first run on a combined 
F606W+F814W image to provide a deep 'master' list of 
objects. It was then run on each of the individual F300W, 
F450W, F606W, and F814W images, the resulting catalogs 
being matched and merged with the master version. In all 
cases, a detection threshold of 1.3(t (where a is the standard 
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deviation of the background noise distribution) and a min- 
imum area after convolution with a 0.3 arcsec top-hat filter 
of 0.05 arcsec'^ (30 connected pixels) were adopted. This re- 
sulted in an average of 200, 420, 510, and 420 objects being 
detected on each WFC chip in F300W, F450W, F606W, and 
F814W, respectively. 

SeXtractor computes several different types of magni- 
tude for each detected object but we worked solely with 
the corrected isophotal (MAGJSOCOR) values as they ap- 
peared to be the most stable at the faintest limits. This par- 
ticular magnitude provides an estimate of the total fiux from 
an object by taking the light measured within the thresh- 
old isophote and then, by assuming a gaussian profile, mak- 
ing a small (< 5%) correction for the light lost outside this 
isophote. As a check, these magnitudes were compared with 
the 'total' magintudes available in catalogs generated at the 
Space Telescope Science Institute using the FOCAS pack- 



age dValdez 19931 ). With the exception of the F300W band 
(see below), excellent agreement was found between the two, 
both in zero point and scale. The scatter observed between 
the two sets of magnitudes was also consistent with the ran- 
dom photometric errors expected from photon statistics. 

It was noticed that SeXtractor and FOCAS both had 
difficulties in reliably detecting and measuring objects on 
the F300W images. A visual inspection revealed that many 
of the detected objects were spurious and yet a consider- 
able number of real objects were missed. This is undoubt- 
edly due to the problems experienced in fiat-fielding the 
F300W data and removing scattered light. A conspicuous 
checker-board pattern peculiar to the drizzling scheme was 
also present around the edges of the frames. These prob- 
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lems were dealt with by using the IRAF APPHOT routine 
to measure aperture magnitudes at each of the positions in 
our master catalogue, setting the aperture diameter equal to 
the major axis length computed for each image on the com- 
bined F606W+F814W frame. It was these F300W aperture 
magnitudes that were used in the selection process. 



3 RESULTS 



Figure 2. Histogram of Absolute Magnitudes. These become 1.4 
magnitudes brighter for a qo = 0.05 universe. 



Using the above colour and magnitude criteria, we searched 
our SExtractor catalogue for candidate high redshift objects. 
We found 9 candidates at F450W<25.7 which satisfied our 
colour criteria. Visual inspection of all these candidates led 
us to reject one of these because it lay near the edge of one of 
the CCD chips and was severely contaminated by spurious 
flux at the edge of the dithering pattern. 

Basic data on the objects is presented in Table 1. Images 
of the objects are shown in Figure la, using the combined 
F606W-I-F814W frame. Corresponding F300W images are 
shown in Figure lb, demonstrating the lack of any detec- 
tions, and thus the presumed presence of the Ly break in 
this band. 

We also used a similar search technique to find can- 
didate objects at still higher redshift by looking for a flat 
F606W -F814W colour combined with F450W-F606W > 2 
and little or no flux in the F300W band. This would select for 
objects with 3.9<z<5, though with additional complications 
from the Lya forest absorption which becomes stronger. No 
such objects were found. 



Figure 3. Histogram of Object Sizes. These increase by a factor 
of 2 for a qo = 0.05 universe. 



4 THE PROPERTIES OF CANDIDATE Z~3 
GALAXIES 

We found 8 objects on the 3 WFC chips likely to be in the 
range 2.6<z<3.9. We now calculate the number-density, ab- 
solute magnitude and physical size of these objects; in doing 
so, we assume them all to be at the mid-point of this red- 
shift range (ie. z = 3.25) and adopt an Ho=100, qo=0.5 
(i.e. Einstein-de-Sitter) cosmology except where stated oth- 
erwise. 

A total of 8 objects within the 4.7 sq. arcmin area of 
the HDF gives a surface density for these objects of 1.7 
arcmin"^. At matching magnitudes we find a similar sur- 
face density of 0.4 arcmin"^ to that of the R~ 25 3<z< 3.5 
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Figure 1. Images of the Objects. The objects are shown in order from left to right then down the page, starting with 2P1 in the top 
left and ending with 4P2 in the bottom middle. All images are 4"square. a Coadded F606W + F814W image, b F300W image showing 
lack of flux in this band for all objects. 
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galaxies discussed by Steidel et al (steidel et al. 1996). For 



the assumed redshift and cosmology, the number density of 
high redshift galaxies in the HDF corresponds to a comov- 
ing number density of 0.005 h^Mpc"^ (where h= Ho/100); in 
comparison, nearby bright (L > L*) galaxies have a number 



density of 0.015 Mpc"^ ( |Parkes et al 1994| ). The similarity 
of these two numbers suggests that we may finally be seeing 
the progenitors of nearby bright galaxies. The corresponding 
number for a qo=0.05 universe is 0.0009 h''Mpc^'^. 

Absolute 'blue' magnitudes were determined from the 
apparent F450W magnitudes using the K-correction cal- 



culated by Cowie et al. (1994) for a star-forming ('spi- 
ral') galaxy at this redshift. This is based on local galaxy 
SEDs from Coleman, Wu & Weedman (1980) in the opti- 
cal, and Mobasher, Ellis & Sharpies ( 1986| ) in the infrared. 
We also apply the -1-0.4 magnitude correction between the 
F450W and standard B band appropriate for a flat spec- 
trum object using the information provided by Holtzman et 



al ( 1995 ). These absolute magnitudes must be treated with 
great caution as neither the applicability of this particular 
K-correction to our sample of high-redshift candidates nor 
the required knowledge of the far-UV spectral energy dis- 
tribution of different present-day types are well determined. 
Indeed the K correction is similarly large and uncertain for 
all the filters discussed here since the objects lie at such a 
high redshift. This notwithstanding, it is of much interest 
that the values found are at the upper end of the luminos- 
ity function of local galaxiesQ (see Figure 2). In Figure 3 

t We note that this calculation is very sensitive to qo with our 
values being about 1.4 magnitudes brighter if a value of qo=0.05 
is adopted. 



we show the sizes of the major axis of the objects measured 
at an isophote 3 magnitudes fainter than the peak surface 
brightness. In contrast to the luminosities, these appear to 
be significantly lower (median diameter ~4 kpc) than those 
of their present-day counterparts. A similarly small-sized 



population of sources was noted by Dressier et al. ( 1994 ) 
in a galaxy cluster found serendipitously on a WFPC-2 im- 
age and thought to be associated with a QSO at z=2.2. 

The object selection is sufficiently bright in the longer 
wavelength filters that the images are deep enough to allow 
examination of the morphologies of the selected galaxies. 
Inspection of the images in Figure la immediately demon- 
strates that all of these objects are disturbed systems, with 
multiple nuclei, jets, plumes, and asymmetric fiux distribu- 
tions. Indeed, the structures observed are similar to inter- 
acting/merging galaxies observed locally e.g. Clements et al. 



( 19961 ). It has been known for some time that interactions 
and mergers between galaxies can trigger bursts of star for- 



mation (Joseph & Wright 1985). These morphological signa- 
tures are also common at intermediate redshifts where HST 
and other high-resolution observations have revealed objects 
of this type amongst the blue 'Butcher-Oemler' population 



in clusters (Couch et al. (199-4), Dressier et al. (1994)), the 



faint blue field population (Glazebrook et al 1995b), and 



mjy radio sources (Windhorst et al. 1995) 



Theoretical models of galaxy formation and evolution 



(e.g. Kauffman (1995)) currently favour the idea that the 
galaxies we see today were formed as a result of the merg- 
ing/accretion of a number of subunits. If this is the case then 
each such accretion event is likely to be accompanied by the 
vigorous star formation we see in local mergers. The dis- 
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Name 


RA(J2000) 


Dec(J2000) 


F450W 


F606W 


F814W 


Abs. B Mag 


Size (kpc) 


2P1 


12:36:45.29 


62:13:47.9 


25.6 


25.6 


26.1 


-21.6 


4.5 


2P2 


12:36:53.34 


62:13:30.4 


25.3 


25.1 


25.4 


-21.9 


4.2 


2P3 


12:36:51.97 


62:14:34.8 


25.6 


26.1 


26.6 


-21.6 


4.0 


2P4 


12:36:55.07 


62:13:49.1 


25.4 


25.7 


26.2 


-21.8 


3.0 


3P1 


12:36:49.85 


62:12:44.6 


25.4 


25.9 


26.5 


-21.8 


5.0 


3P2 


12:36:58.20 


62:12:17.7 


25.5 


25.8 


26.2 


-21.7 


3.9 


4P1 


12:36:45.25 


62:11:54.5 


23.9 


23.8 


24.2 


-23.4 


10.9 


4P2 


12:36:41.58 


62:12:39.9 


24.5 


25.0 


25.4 


-22.7 


3.1 


Table 1. 


Candidate Primeval Galaxies 


in the HDF. 









All magnitudes are on the ST magnitude system. No galaxies were detected in the F300W filter, giving a 2 cr magnitude limit >~27.7. 
Names follow the pattern chip no. P object no. so that the first object detected on chip 2 is 2P1. Absolute mag-nitude s are calculated 



assu ming a K correction appropriate for spiral galaxies and the appropriate magnitude shift from F450W to B given in (Holtzman ct al 



1995 |). y alues are given for Ho = 100 km ^Mpc— 1 and qo=0.5. For qo = 0.05, absolute magnitudes get brighter by 1.4 mags, and 
increase by a factor of 2, while for Ho=50, corresponding figures are 1.5 mags brighter and 2 times bigger. 



turbed, multiple-component nature of the HDF high redshift 
galaxies presented here suggests that this is a wide-spread 
occurrence at these very early epochs. The close match be- 
tween the number densities of these objects and local galax- 
ies might indicate that we are at last seeing the epoch where 
most of the galaxies we see today were being assembled from 
such subunits. Alternatively, if there has been considerable 
number-density evolution in galaxies at lower redshifts (0.5 
- 1), as might be indicated by the faint- blue-galaxy popula- 
tion, then we might just be seeing occasional starbursts in 
a population of more numerous small galaxies. 

The effect of viewing these objects in their restframe 
UV (~ 1670 A) is also another factor that must be consid- 
ered. We might just be seeing the hottest, brightest regions 
of active star formation in a larger, lower surface brightness 
structure that we cannot see. Observations of the HDF at 
IR wavelengths, particularly those conducted with the new 
NICMOS camera soon to be installed on HST, will be in- 
valuable in addressing the problem of the observed restframe 
wavelength. For example, the centre of the K window cor- 
responds, at these redshifts, to ~ 5100 A, conveniently in 
the restframe optical emission. Observations at these wave- 
lengths from the ground will also be useful in better con- 



straining the spectral nature of these objects; we shall be 
conducting these in the forthcoming months. 

5 CONCLUSIONS 

We have presented results on a sample of objects very likely 
to lie at 2.6<z<3.9. We find that these objects have sim- 
ilar co-moving number density to nearby bright galaxies. 
However they are significantly more luminous and smaller 
than nearby galaxies. Almost all of our high redshift objects 
show signs of a disturbed and multi-component morphol- 
ogy, which we interpret as being due to merging processes 
taking place during the assembly of a galaxy. 
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